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The packing problem of long thin filaments that are injected into confined spaces is of fundamental 
interest for physicists and biologists alike. How linear threads pack and coil is well known only 
for the ideal case of rigid containers, though. Here, we force long elastic rods into flexible spatial 
confinement borne by an elastic shell to examine under which conditions recently acquired knowledge 
on wire packing in rigid spheres breaks down. We find that unlike in rigid cavities, friction plays 
a key role by giving rise to the emergence of two distinct packing patterns. At low friction, the 
wire densely coils into an ordered toroidal bundle with semi-ellipsoidal cross-section, while at high 
friction, it packs into a highly disordered, hierarchic structure. These two morphologies are shown to 
be separated by a continuous phase transition. Our findings demonstrate the dramatic impact of 
friction and confinement elasticity on filamentous packing and might drive future research on such 
systems in physics, biology and even medical technology toward including these mutually interacting 
effects. 
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Introduction. — Dense filament packing can be found 
in various natural systems. A well known instance is the 
injection and subsequent coiling of long DNA in glob¬ 
ules and viral capsids (e.g., m\)- A similar technique 
has been harnessed by neurosurgeons for the minimally 
invasive treatment of saccular aneurysms, into which de¬ 
tachable platinum wires are fed to initiate occlusion [4]. 
Extremely dense fiber packing can also be observed for 
example in gland thread cells of hagfish [5] . Albeit their 
morphology has recently been unraveled [6], the under¬ 
standing of morphogenesis—the process of shape transfor¬ 
mation and development—remains vague. In particular, 
the role of macroscopic material parameters is little un¬ 
derstood. 

Several experimental and numerical studies [THU] re¬ 
vealed how long wires form loop patterns and alignment 
between contacting segments when injected into rigid 
two-dimensional containers, depending on friction, plastic 
yield point and the precise insertion setup. In rigid three- 
dimensional cavities HHE], on the other hand, friction 
has been reported to have negligible impact on the pack¬ 
ing process. In a recent study on morphogenesis of elastic 
ring filaments growing in flexible shells m, however, we 
discovered four morphological phases strongly dependent 
on friction, which brought it back into play in 3D. This 
naturally begs the question to which degree the knowledge 
previously acquired on the thread injection and packing 
problem in rigid cavities is applicable to less idealized, 
deformable ones, such as cell walls, vesicles, or arterial 
walls. In this paper, we perform a change of topology 
from ring-like (as in Ref. m) to linear threads, which 
are more relevant in Nature and biomedical applications, 
to answer this question quantitatively using numerical 
simulations and simple table-top experiments. We show 
how confinement elasticity completely alters the packing 
process of a thin wire that is fed in. In strong contrast 


to rigid cavities, friction becomes the key macroscopic 
material property beyond a critical point, determining 
the packing process by giving rise to the spontaneous 
emergence of two different morphologies. We characterize 
them by comparing numerical measurements of their en¬ 
ergetic and structural properties, finding that frictional 
forces lead to a highly disordered, hierarchic, crumpled 
structure governed by power laws, whereas absence of fric¬ 
tion yields a dense toroidal bundle with semi-ellipsoidal 
cross-section. We further show that the transition from 
rigid to flexible cavities as well as the transition from weak 
to strong friction in flexible cavities is continuous and 
accompanied by spontaneous symmetry breaking, and we 
identify the associated order parameter. These results 
challenge researchers in the many fields where filamentous 
packing is relevant to carefully include the effects of fric¬ 
tion and elasticity in their studies, and might trigger new 
efforts for instance in aneurysm coiling research, which 
has previously neglected these. 

Model. — For the theoretical considerations and nu¬ 
merical simulations, we model the thin wire by an ex¬ 
tensible, intrinsically straight and untwisted Kirchhoff 
rod m- Its centerline, a space curve x(s), is parame¬ 
terized by its arclength s G [0,1/]. In its deformed state, 
an orthonormal director frame {di}i=i^ 2,3 specifies the 
cross-sectional orientation along the curve, with the third 
one being the unit tangent ds = ||dsx||. Its rate of 

change d^di = k x d^ defines the Darboux vector k with 
director components ki = k • d^. Assuming a homoge¬ 
neous, isotropic, linearly elastic material with Young’s 
modulus and Poisson’s ratio the elastic potential 
energy of the wire reads 

C/w = f / + GJkl + ds, (1) 

^ 7o 
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in which k = ^Jk\^k 2 denotes the curvature, k^ is the 
twist per unit length, whereas e = ||5sx|| — 1 is the axial 
Cauchy strain due to compression or tension. For an 
invariant circular cross-section with radius r, A = irr^ is 
the cross-section area, / = irr^/A the second moment of 
inertia, and J = 2/ the polar moment of inertia. The 
shear modulus reads G = E^/2{1 + z^w)- 

Analogously, a Kirchhoff-Love shell with middle sur¬ 
face ft and uniform thickness t is used to represent the 
flexible hull. In linear elasticity (Young’s modulus and 
Poisson’s ratio z^s), the elastic potential energy of the shell 
reads [18] 

= l f {Maijaki + BpijM dn, (2) 

^ Jq 

(3) 

where M = E^t/{1 — z^g) the membrane rigidity and 
B = Mt^ j VI the bending rigidity. Einstein summation 
over repeated indices i^j^k^l G {1,2} is used, with sub¬ 
scripts (superscripts) denoting covariant (contravariant) 
coefficients. The membrane strains aij = (a^j — aij)/2 
and the bending strains ^ij = hij — bij are the covariant 
coefficients of the change of first (a) and second (6) fun¬ 
damental forms from the initial stress-free configuration 
(barred symbols) to the current one (bare symbols). 

For a realistic treatment of body contacts between any 
combination of types (wire-wire, wire-shell, shell-shell), we 
use a contact model with dry stick-slip friction analogous 
to Refs. [T2j[T4l[T6| in the simulations. Repelling Hertzian 
contact forces are exchanged when two body volumes 
overlap, whereas the tangential contact forces are set to 
obey Coulomb’s law with static and dynamic friction 
coefficients /ig and /id, respectively. 

We minimize the elastic energies U = 1/w + with the 
finite element method and integrate Newton’s equations 
of motion with a common predictor-corrector scheme of 
second order with adaptive time-stepping. Subcritical 
viscous damping forces are added to allow the system to 
stay close to static equilibrium at all times. Details on 
these numerical models and their implementation can be 
found in Refs. naEoj. 

An intrinsically straight thread is radially injected into a 
thin shell whose unstrained equilibrium configuration is a 
sphere with radius R with a small opening with diameter 
2r through which the thread is fed in. The confining 
shell is held in place by imposing zero displacement on 
a narrow rim about the entrance hole. A small random 
transverse deflection is initially added to the wire to 
break the rotational symmetry about the feeding axis. 
In all simulations, the following parameters are fixed: 
r = 0.5 mm, E^ = 1 GPa, z^w = = 1/3. To reduce 

the influence of friction on the phase space to a single 
dimension, we always use a fixed ratio jUd/Vs = 0-9. The 
wire is pushed in at a constant speed of 0.5 m/s, which is 



FIG. 1. Packing evolution in flexible confinement. Series of 
simulation snapshots at /Xg = 0 (top row) producing aligned 
coils, and at /Xg = 1/2 (bottom row) producing disordered 
crumpled packings. The remaining parameters are /? = 20, 
^ = 200, £ = 100. All images to scale. 


slow enough for inertial effects to have negligible impact 
on the outcome. 

In summary, the system is effectively characterized by 
the four dimensionless and independent control parame¬ 
ters p = R/r, ^ = R/t, 5 = E^/Es and /Xg if it is close to 
static equilibrium—a parameterization which is equiva¬ 
lent to that introduced in Ref. m- The slenderness ratio 
^ and the cavity flexibility 5 generalize the rigid spheres 
considered in preceding related studies pISHTS] . 

Packing evolution and morphologies. — The typi¬ 
cal packing evolution obtained in sufficiently elastic or thin 
shells is displayed in Fig. for increasing thread lengths 
I = L/R, and a corresponding movie is provided in the 
supplementary material (see video_sl.mpg). In absence of 
friction, a highly ordered packing pattern is observed in 
which the wire bundles into a tight toroidal coil (compara¬ 
ble, e.^., to how microtubules bundle up in erythrocytes 
[21], or to DNA spools in bacteriophage heads [T]) that 
continues to grow and stretch the confining shell as more 
thread is injected. Eventually, both upper and lower faces 
of the shell are fully flattened, turning it into the con¬ 
vex hull of the enclosed coil. Similar forms have been 
experimentally obtained by enclosing elastic nanotubes 
and nanowires with emulsion droplets and polymer shells, 
which were then intentionally contracted [22] [23| . Monte 
Carlo simulations at finite temperature (2ll[2l] have like¬ 
wise indicated that soft vesicles deform into obloids when 
enclosing a fluctuating polymer chain whose persistence 
length grows much larger than the vesicle diameter. 

When friction is activated, however, the situation is dra¬ 
matically changed: Tangential sliding is hampered, which 
lets the wire tip poke the surrounding shell significantly. 
The compressive forces acting on the inserted thread are 
much higher in consequence, and soon let it buckle out 
of the coiling plane to form a more three-dimensional 
packing process with frequent spontaneous loop reorienta¬ 
tions, yielding a crumpled structure that tends toward a 
spherical globule at high packing density. Due to the thin 
shell’s tensile flexibility, frictional forces let the pushing 
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FIG. 2. Packing morphologies in flexible confinement, (a) 
Toroidal coiling at weak friction (/Xg = 0). (b) Crumpling at 
strong friction (/Xg = 1/2). The further simulation parameters 
are p = 20, ^ = 10, s = 10^, and the packed wire length 
is I = 800. In color, the dimensionless curvature is shown. 
Images to scale, except that the wire radius is halved to reveal 
the inner structure. Movies for (a) and (b) are available in 
the supplementary material (video_s2.mpg and video_s3.mpg). 
(c-e) Experimental realizations of the two morphologies in 
transparent and opaque rubber balloons, (c) r = 1 mm, 
R— 23 mm, I ^ 950. (d) r — 0.5 mm, R — 23 mm, /Xg ~ 0.55, 
I ^ 650. (e) r = 0.5 mm, R = 18 mm, /Xg 1.2, I ^ 750. (f) 
Cross-section of (a) with superposed ellipses and geometric 
parameters. 


wire locally drag its confinement along a small distance, 
making static friction much more relevant in deformable 
than in rigid cavities. 

For validation and demonstration, we have reproduced 
these two morphologies in table-top experiments, three 
examples of which are shown in Fig. c-e). Polycaprolac¬ 
tam wires (r = 0.5,1 mm, up to L = 22 m) were manually 
pushed through a straight steel pipe into customary in¬ 
flatable balloons made of natural rubber {R= 18, 23 mm, 
t = 0.25 mm). The friction coefficient was determined 
by placing pieces of wire on pieces of the rubber shell 
and measuring the tilt angle at which the wire started 
to slide. To reduce friction, we coated the wires with 
an acrylic dye and additionally used a silicone lubricant. 
These experiments yielded the same two morphologies as 


our numerical simulations. 


Analytical model for the coiled morphology. — 

To quantitatively characterize the coiled morphology 
at low friction, we measured the geometrical properties 
of cross-sections, of which a typical example with fully 
stretched lateral faces of the thin shell is displayed in 
Fig. Ha ,f). Within the reasonable parameter range ex¬ 
amined here (20 < p < 40, 10 < ^ < 200, 1 < e < 10^), 
if the shell is flexible enough, the cross-sections of the 
toroidal bundle are very well approximated by two half 
ellipses, with major toroidal radius i7t, and two minor 
ellipsoidal radii Rx^ Ry^ as labeled in Fig. j^f). While 
the evolution of these three radii with the packed thread 
length depends on the system parameters p, 5, the com¬ 
mon functional form is well approximated by a power-law 
scaling Rt/R ^ RxjR ^ RyfR ^ P in all ob¬ 
served cases within a certain range of the rescaled thread 
length I = L/R. This is demonstrated in Fig. Ha) on 
an example simulation at p = 20, ^ = 10, 5 = 10^ for 
which a = 0.04 ± 0.01, P = 0.35 ± 0.01, 7 = 0.46 ± 0.01, 
indicating that while the minor radii grow quickly, the ma¬ 
jor radius Rt remains approximately constant, a = 0 is 
indeed expected in the flexible cavity limit M ^ Egt ^ 0 
(x.e., ^ 00), since 1/Rt is the maximum curvature 

in the coil, the square of which is minimized according 
to the principle of minimum energy. This high degree 
of order allows for an analytical approximation following 
the ideas of Purohit et al. [26]. Their geometrical model, 
which was originally devised for the description of DNA 
coils in viral capsids, has also been successfully applied to 
ordered wire packing in rigid cavities [HUS]. It is based 
on the approximation that thin threads coil such that 
their binormal vectors are always parallel to the main coil¬ 
ing axis X, and that their radius of curvature about this 
axis is maximal, resulting in an empty cylindrical region 
as it is observed here. Assuming that the centerlines of 
individual strands in the coil are separated by a distance 
d{L) > 2r on average, the number of windings along the 
X axis is given by w{y) = 2Rx^Jl - {[y - Rt]/RyY/d{L) 
for y > Rt. The packed thread length follows as 


L = 


2 


L 


Rt-\-R 


y 

2TTyw{y) dy = 


2TrcRxRy 


(4) 


where c = 4/3 + 7rp and p = Rt/Ry. This fixes the 
average strand separation d{L) given that L, R^, R^, Ry 
are known. Analogously, the wire bending energy Uh = 
1/2 E^Ik^ ds reads 


Ub{L) = 


V^d{L) 


! 

JRt 


Rt+Ry I 


R’iy) dy. (5) 


Using Eq. the unknown d{L) can be eliminated to 
yield 


Ub{L) = 


2LEy,I 

cR? 

uriy 


q In 


1 + y 


|y-i 


(6) 
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where q = — p^. Note that Eq. ^ is real-valued even 

for p > 1. U\) is monotonically decreasing in p, and thus 
the wire favors toroidal configurations with large major 
radius Rt but small minor radii Rx^ Ry. This tendency of 
the wire competes with the shell deformations necessary 
to adopt such a shape. 

Quantitative comparison of the two morpholo¬ 
gies. — In Fig. [^c), the dimensionless bending ener¬ 
gies of both morphologies are compared. Since Uh ^ 
p^E^ljR in the rigid shell limit [14], we rescaled it to 
Uh •= U\yr‘^/ REy^I . The numerical measurement of the 
coiled morphology is very well approximated by Eq. (§, 
with only a slight overestimation at early packing. The 
crumpled morphology simulated at /ig = 1/2, on the 
other hand, emerges by elastic bifurcation, which can 
be recognized in Fig. |^c) by the early discontinuity in 
U\y. It exhibits a clear power-law trend I/b ^ with 
S = 1.192 ± 0.006, which hints at a hierarchic inner pack¬ 
ing structure: The more wire is injected, the more spatial 
freedom is limited, resulting in higher bending curvature 
in the newly formed loops. The bulges visible in Fig. |^e) 
provide some intuition on this phenomenon. Such a power- 
law scaling of the bending energy discriminates the dis¬ 
ordered phase in elastic shells reported here from that 
obtained in rigid spheres with intrinsically curved wires 

[El* 

To gain further insight into the packed structure, we 
consider the total curvature K = k ds, which reads 


K(L) 


2 f^t+Ry 

7m L 


( 7 ) 


in our coiling model (recall c = c{L) as defined below 
Eq. @). This excellently approximates the numerical 
measurement with a slight initial overestimation, as shown 
in Fig. id). In the hierarchic crumpled morphology 
at /is = 1/2, we find a power law K ^ with A = 
1.083 ± 0.004 consistently with the energetics. Again, 
a pronounced discontinuity in K marks the beginning 
of its existence at the bifurcation point. The hierarchic 
nature of the crumpled structure is manifest also in the 
number of contacts between discrete wire segments {in 
silico: finite elements), as Fig. [^b) shows. After the 
two morphological phases have bifurcated, it follows a 
power law N E with r = 1.40 ± 0.01 at /ig = 1/2. On 
the contrary, low friction leads to dense alignment with 
hexagonal packing in the coiled morphology, hence N ^ I 
once a few coil windings are established. 

Hierarchic disorder. — Above, we have termed the 
disordered morphology hierarchic^ based on the observed 
power-law scaling of the contact number, bending energy 
and total curvature. Now we consider the geometric 
structure of the packed wire in additional detail to put 
this statement on more quantitative ground. 



FIG. 3. Quantitative comparison of the two morphologies, (a) 
Scaling of the shape of the coil with the inserted wire length. 
Data obtained by fitting semi-ellipses to the cross-sections as 
shown in Fig. [^f). (b) Number of discrete wire self-contacts 
N. (c) Non-dimensionalized bending energy Uh- (d) Total 
curvature K. Symbol and line legends in subplots (c) and 
(d) hold for (b-d). Data from simulations at /? = 20, ^ = 10, 
£= lOT 


In rigid spheres, the distribution of bending energies in 
disordered wire packings is well fit by a lognormal proba¬ 
bility density, with the exception of a larger tail toward 
low values na. The same holds here for the crumpled 
morphology in elastic shells, albeit with a different origin. 
In rigid spheres, such hierarchic disorder needed to be 
introduced by artificially pre-curving the wire, whereas 
tight membrane confinement induces this hierarchy natu¬ 
rally even with straight threads. Figure [^a) shows that 
the logarithm of the normalized local squared curvature 
follows a normal distribution: 

with probability density function 


/(^) = 


1 


exp 


1 f X — 

2\^ j 


(9) 


Our numerical data for the crumpled morphology is best 
fit by Eq. ^ with a mean value of T = —0.29 and a 
standard deviation of cr = 0.83. 

Lognormal distributions are typical for hierarchic events. 
In crumpled thin sheets, for instance, the ridge lengths 
are lognormally distributed, see Ref. m and the liter¬ 
ature cited therein. For the present crumpled wires in 
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FIG. 4. Hierarchy in the disordered morphology, (a) Distri¬ 
bution of the squared curvature j (b) Distribution of 

loop sizes A/r. Asterisks denote numerical data points, solid 
lines are best-fitting lognormal probability densities. 

elastic shells, the relevant local geometrical features are 
individual loops formed in the packing process. A wire 
loop is easily defined in two dimensions [7], but in 3D it 
is somewhat less obvious. We define it here as the wire 
segment between two neighboring local maxima of the 
squared curvature along the wire centerline, and denote 
its length by A. As shown in Fig. |^b), the distribution 
of loop sizes in the crumpled morphology is lognormal in 
good approximation: In (A/r) ^ A/'(x, cr^) with x = 3.0, 
cr = 0.48. 

Order-disorder transition. — To complete the quan¬ 
titative characterization of the two morphological phases, 
an order parameter is required that discriminates them rig¬ 
orously. What defines the order of the coil at low friction 
is alignment between the individual loops. The hierarchic 
packing pattern in the crumpled morphology emerges due 
to frequent three-dimensional loop reorientations, which 
break this alignment. Denote by n the unit vector point¬ 
ing in direction of the wire’s principal axis of minimal 
moment of inertia. The curvature of the wire about n 
is given by the triple product = n • (d^x x Q^x), 

whose average sign S = 1/L sgn ds measures 
the fraction of the wire turning in either direction about 
that axis. S vanishes iff left- and right-turning wire seg¬ 
ments are balanced {e.g.^ when the loops are isotropically 
distributed, degenerate n), whereas it takes one of its 
extreme values ±1 iff the coil never changes orientation. 
Since this coiling direction is initially selected at random, 
the sample average of S is zero. The decisive non-trivial 
quantity is thus |5'|, and = 1 — |5'| is an order param¬ 
eter that discriminates ordered coiling from disordered 
crumpling. A closely related quantity has already served 
to distinguish morphologies of growing elastic rings con¬ 
fined in spheres m- As shown in Fig. D vanishes in 
the coiled phase observed in relatively stiff confinements 
even with very strong friction, and in flexible shells if 
the friction coefficient is sufficiently small. Evidently, 
the crumpled phase {D > 0) requires a flexible shell and 
moderate friction to emerge, and the phase transition is 



FIG. 5. Order parameter D for the coiling-to-crumpling tran¬ 
sition. (a) As a function of the conhnement flexibility s at 
p — 20, = 30, ps — 1. (b) As a function of the friction 
coefficient /is at p = 20, ^ = 10, £ = 10^. The inset shows a 
magnihcation about the critical point. The solid gray lines 
are used as a guide to the eye. Error bars indicate standard 
errors from 10 independent simulations. Data maximized over 
I G [5,50]. 


continuous. 

Conclusion and significance. — Our findings accen¬ 
tuate the fundamental difference between rigid and flex¬ 
ible cavities for the wire packing problem. Friction is 
responsible for a morphological phase transition from or¬ 
dered coiling to a disordered hierarchic structure. In rigid 
spheres, such disorder needed to be introduced by artifi¬ 
cially pre-curving the wire [14], whereas in the tight flexi¬ 
ble confinement considered here it occurs spontaneously 
even with straight rods. This insight might contribute 
to explaining the high degree of order observed in hag- 
fish gland cell thread morphogenesis. It has far-reaching 
implications for biomedical applications such as the sur¬ 
gical occlusion of cerebral aneurysms, where the friction 
coefficient is estimated to be well above the critical value 
reported here and where common simulations ignore de¬ 
formations of the arterial walls (e.p., [28]). While friction 
can cause problems during injection [29], our results sug¬ 
gest that it is nevertheless highly desirable in order to 
obtain a broader spatial distribution of loops, which is 
known to assist in stabilization [28] . 

Elastic threads with circular topology develop four dif¬ 
ferent morphologies when flexibly confined m- We have 
found in this letter that only two prevail when linear 
threads are injected instead. The reason is that elastic 
rings possess an instability when confined in spherical 
shells that linear threads do not possess: They are (ini¬ 
tially) in full contact with the inner surface of the shell 
and can either buckle away from it or not. They remain 
on the two-dimensional surface in rigid spheres and are 
enabled to pack in a three-dimensional fashion in flexible 
shells. Linear wires, on the other hand, are never in full 
contact with the shell, due to their moment-free ends and 
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the injection at a right angle. They hence avoid being 
trapped in the energetically unfavorable surface packing 
configuration, avoiding the instability in the first place. 

An open question is the sensitivity of the packing behav¬ 
ior to material nonlinearity. It could also be worthwhile 
to determine the functional dependence of the exponents 
and critical points reported here on the system parameters 
to obtain a full image of the four-dimensional parameter 
space. 
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